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1. Introduction 

In this paper, we are mainly concerned with the mathematical analysis of the evaporation of a single 

9 

droplet in a gas, in the continuation of our previous work—'. 

Experimental studies around this sub ject are of course important for industrial purposes. Let us 
refer for instance to the work a'^ l ' l " l I "-^ l 

One of our main interest in this paper is to analyze the time evolution of the droplet radius, a study 
that we began in our previous worli^. Let u s just mention here that the experimental evolution of 
this radius is well known as the d where d denotes the diameter of the droplet, see below 

for more details. 

Let us recall the standard physical framework for this evolution. 

The evaporation of a single droplet in a gas involves simultaneous heat and mass transfer processes. 
In particular, heat from evaporation is transferred to the droplet boundary by conduction and 
convection, while vapor is carried by convection and diffusion back into the gas stream. Evaporation 
rate depends on the pressure, temperature and physical properties of the gas, the temperature, 
volatility and diameter of the drop in the spray. 
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To fix the ideas, in the experimental study^ of a single droplet evaporation performed by the 

LCSR (Combustion Laboratory of the University of campus from Orleans, France), the droplet is 

suspended from a silicate tube. The elliptic shape of the droplet is assimilated to a sphere of equal 

volume. Important quantities of interest for these experiments are time evolution of the droplet 

radius, as well as classical quantities such as mass fractions or temperatures of the liquid and gas. 

In the experimental studies performed above, the so-called d? law is used to simplify two-phase 

fluid models and then propose adequate numerical schemes. This law simply states that the time 

d"^ ■ ■ 

evolution of the radius behaves as — in time flow, and is purely phenomenological. 

"o 

Our purpose in this paper is exactly in the opposite sense. We start from phenomenological fluid 
(mixtures) PDE modeling the drop evaporation process, compute the time-evolution of the drop 
radius, and then deduce other quantities of interest such as mass fractions of the liquid and gas. In 
particular, our numerical experiments are in good agreement with this phenomenological d"^ law, at 
least for small time evolution. 

Our framework is therefore as follows: we consider a droplet initially represented as a single compo- 
nent mixture (liquid chemical specie 1) while the surrounding gas at time i = is made of only one 
(gas) chemical specie, say 2. 

During the evaporation process, the liquid vapor is transferred into the gas, while by condensation 
at the droplet surface and then by diffusion, gas chemical specie 2 appears inside the droplet. 
We make the important simplification that the moving interface between the droplet and the sur- 
rounding gas (i.e. between the two species) is spherical, with radius R = R{t) evolving in time. 
Let pG (resp. pl) denote gas density (resp. liquid density), and vg (resp. vl) denote the gas velocity 
(resp. liquid velocity). Then, one has the classical overall continuity and momentum conservation 
laws 

dtPk + div(pj.Wfe) = (1.1) 

PkdtVk + PkVk-^Vk = -Vp. (1.2) 

Above, subscript k refers to the gas G or to the liquid L, depending on whether one considers the 
gas or liquid.p the state equation of the gas. 

Let Yli,Yl2 (resp. ^01,^02) the mass fractions of the hquid (resp. gas) obtained after diffusion of 
species in the surrounding gas. Therefore for two species, one has 

Ygi+Yg2=Yli+Yl2 = 1. 

Along with equation (jl.ip . we have to add the equation giving species conservation. So for the liquid, 
we have 

PLdtYLk + PLVL-^YLk + ^HpL^LkVLk) - -pL/(i^Lfe), A: = 1, 2, (1.3) 

Y^k denoting mass fraction of the liquid, and / being a continuous function modeling a friction or 
a resistance for the drop. 

We assume that the liquid speed is so small that is can be settled to 0. Equations (jl.ip and (|1.3p 
can then be written under conservative form as 

dt{p~g) + diY{p~gv)^F{~g) (1.4) 
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or 



d,u+-^{f{u)) = F{u), (1.5) 

in a system of particular coordinates. 

If r is a curve of discontinuity of u, then one has 

[/(«)] = Nf, (1.6) 

where [.] denotes the jump of the inner quantity, s — ^ is the speed of discontinuity along F. The 
jump relation p.6p is known as Rankine-Hugoniot condition. It merely means that discontinuities 
cannot be completely arbitrary. The above considerations are all classical facts^ 
In the case of our droplet, in order to find interface condition at the droplet surface, i.e. for r — R(t), 
it is sufficient to use ()1.5|) and ()1.4p in polar coordinates, getting 

[pgv] = [pg]^. (1.7) 

Thus taking 5=1, one has 

dR 

[PG - Pl)-j^ = Pgvg - Plvl, 



that is also with vl — 



dR\ dR 



Taking g = y in (|1.7|) . Y denoting the mass fraction of the liquid or the gas after diffusion, we get 

dR 

{pGYck - PLYLk)-rr = PGYckivG + VGk) - PhYLk^Lk 

dt 

and this is equivalent to 

PGYGk{vG - R') + PGYokVGk = -phYLkR' + PhYLkVLk- (1-9) 

Above VQk (resp.Wife) is the speed of sp ecie Gfc (resp.Lfc), k — 1,2. 
Combining relation (II. 9p with Pick's that is 

YgiVgi -Di2VYg1- ^02^02 = --D21VYG2, 

D12 and D21 being diffusion coefhcients, and with equations relating the thermodynamic state at 
the interface r = R{t) 

Yck^KkYLk, fc= 1,2, 
we obtain, for the mass fraction of the liquid Yli, the boundary condition 

R'(t)(Ki - 1) , , , , 

+ K2PGmu)-K/ ^^ = (i-^^) 
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using polar coordinates. 

In our previous worli^, we have made huge mathematical and physical simplifications taking the 
state equation of the gas p as constant in (|1.2p , and considering gas velocity vg as a given function 
of the time t. Thus in our previous work, system (|1.2p was reduced to equation (jl.ip with a 

given vcit). 

In the present work, we consider the full hyperbolic system (|1.2p with an auxiliary state 

equation for the gas given by pi = . This of course extends our previous work, but considering 
such pressure laws has the advantage that we have been able to perform numerical comparisons. 
More general state laws will be studied in a future work. 

Once paif, t) and voir, t) determined, radius R{t) of the drop suspended in the gas will be computed 
through the ordinary differential equation (jl.Sp . Then we shall determine the mass fraction Y^i of 
liquid after evaporation process, through the PDE (|1.3p along with boundary condition (|1.10|) . for 
a given suitable function /. 

For this last purpose, within the framework of weighted Sobolev spaces on initial data and for some 
continuous function / subject to increasing condition, we shall provide an unique local solution for 
the mass fraction Y^i of the liquid. In addition, wc shall show that if the initial condition is bounded, 
then so is our solution. 

In the numerical applications (2nd example) we have chosen the experimental conditions made by 
the LCSR in the study of single drop evaporation and in this case the study of the radius of the 
drop shows us that the graphic associated to our mathematical model presents the same features as 
in the experimental curves. 

Plan of the paper: In Section 2, by using Riemann invariants, we determine the droplet radius. 
This enables to get, in Section 3, the liquid mass fraction, using a variational method. Finally, we 
have presented some numerical simulations in the last Section, which shows that our model is is 
good agreement with experimental simulations, at least for short time. 



2. Hyperbolic system and droplet radius 

The gas velocity voir, t) and its density poir^ t) satisfy the following system, using polar coordinates 

dpG , 1 d 

yr PG'UG) = u 

(2.1) 



ot r'' or 



Setting p(r,t) = r'^pG{r,t) and v{r,t) — VG{r,t), we have 

dp dp dv 

dt dr dr /o r,\ 

dv dv _ I dpi ^^-^i 

dt dr p dr ^ 

where pi{r, t) is an auxiliary function connected to the state equation of the gas p(r, t), by 

dpi ^ 
dr dr 

In (2.2), according to the discussion in the Introduction, we choose the auxiliary function pi{r,t) as 
Pi=pT, 7 > 1. 
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With this choice, we get the foUowing system 



dp dp dv 

ot or or (2 3) 

dt dr dr 



We note that (|2.3p is equivalent to matrix form 

dt\v J dr \v 

A being the (2, 2) matrix 



'^'^ ■ A#f^).0, (2.4) 



V p 
^p''~'^ v 



Eigenvalues of A (characteristics speeds) are given \ = v — c and /z = w + c, c = \/ ^p^~^ . Since 
X < H, system (|2.4p is therefore hyperbolic. Thus there exists two functions W{p,v) and Z{p,v) 
(Riemann invariants) such that 

W(p,v) — constant on — ; — — A, (2.5) 
dt 

dX'^ 

Z(p^ v) — constant on = /i. (2-6) 

dt 

W{p, v) is determined by the system = vector Ri = (1, c/ p) being and eigenvector associated 
to the eigenvalue p. Thus 

W{p,v)=v-^^. (2.7) 
7-1 

Similarly, the Riemann invariant Z(p,v) corresponding to A is given by 

Z{p,v) = v + ^. (2.8) 
7-1 

Functions W{p, v) — W{t, r) and Z{p, v) = Z{t, r) satisfy the following system, equivalent to system 
(EH) 

dW dW 

+ XiW,Z)—^0 
dt dr (o q\ 

dZ , ,dZ ^^-^1 
-+M(M^,^)^=0, 

where A(VF, Z) and /i(VF, Z^ are given by 

A = -(2^)Z+(l±i)W- 

as follows from ((2J)) and ((2?8)) . 

It is well known that a sufficient condition in order that (j2.9p is authentically nonlinear is that 

> and > 0, which is the case here according to (|2.10p . 
Integration along the characteristics defined by 
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— = XiW,Z), x\o) = p 

gives 

4o/3)W=/3+ f\iW{s,X'is)),Z{s,X\s)))ds. (2.11) 
Jo 

Therefore, the solution of the initial value problem 

Wt + X(W, Z)Wr ^ 0, W{0,r) ^Wo{r) (2.12) 

can be written as 

W{t,r) = W^o(^(V/3)(0)) = WoiP), (2.13) 

where p = r - J* A(VF(s, X^{s)), Z{s, X^{s)))ds. 
Similarly we have 

Z{t,r) = Zo[x^o,c.)iO)) = ^o(a), (2.14) 

where 



^(Oa)W=« + / KWis,X'{s)),Zis,X'{s)))ds. 
Jo 

The above considerations lead to the following 

Proposition 2.1. Assume that WQ{f3) < or Z'^la) < 0. Then solution of system (2.3) is defined 
on a finite interval [0, T[. 

□ 

Proof: Differentiation of (|2.5p and (|2.1ip with respect to /3 gives 

= Xf^iW, Z) with X^pit = 0) = 1, (2.15) 

and in the same way 

dX^ 



- ^i^{W, Z) with Xi{t = 0) = 1. (2.16) 

Since A/3 = Xy/Wp + A^Z^ = A^Wo(/3) and = fizZQ{a), from (|2.15p and (|2.16p . integrating 
w.r.t. t along the characteristics yields 

X^it) = 1 + ^ XwW(,iP)dt = 1 + (2.17) 

Xl{t) = l + J^ ^^zZ',{a)dt = l+(^^^Zl,{a)t. (2.18) 

From ((2l7l) . it follows that X}Jti) = for ti = 7 > 0. Similarly X^(t2) = for 

(7 + l)Wo(/3) 

''^hTT)l^>°- 
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Hence ^{t,r) becomes infinite for T = inf {<i,t2}, since ^ = W^/3-|f = 

□ 

From Proposition it follows that 

Proposition 2.2. System h2. 9\) admits an unique solution on [0, T[, for all r G and for 
initial data pc{0,r) — poir) and tiG(0,r) — vo{r) belonging to C^(M+). 

□ 

Concerning the droplet radius, it follows from (|1.8p . that we have the following ode for this radius 

dR{t) _ VG{t,R{t))pG{t.R{t)) 



dt pcit, R{t)) - PL 

We immediately deduce 



i?(0) = Rq. (2.19) 



Proposition 2.3. The Cauchy problem (2.19) has an unique solution R(t) on a maximal time 
interval [0, r*[ with T* < T, given initial data po{r) and vo{r) such that W(5(r) < or ^o(^) ^ 0- 

□ 



3. Liquid Mass Fraction 

The liquid mass fraction Y^i satisfies the conservation equation of specie (|1.3p . which can be rewritten 



as 

^.YL^-^§-^(r^§-n^)+f{YL,) = 0. (3.1) 

We have used polar coordinates, and taken the diffusion constant D12 as being equal to 1. Of course, 
p.ip is equivalent to 

dtYLi - AYli - --^Yli + f{YLi) = 0, for < r < sit), (3.2) 
r or 

where s{t) ~ R{t) denotes the droplet radius determined in section 2. 

The boundary condition at the surface s(t) is given by the Rankine-Hugoniot condition connected 
to the thermodynamic equilibrium, i.e. formula p.lOp . 

Performing the change of variable r = R{t)x, function YLi{t,r) turns to function YLi{t, R{t)x) = 
u{t,x), which satisfies the following initial boundary value (i.b.v.) problem 

dtu- a[t)( l^u+-d^u\ -x?^da:U + fiu) = < a; < 1, i > 0, (3.3) 
V X J R{t) 

\\mix^a+xux{t,x)\ < 00, U2,(t, 1) + fc(t)u(t, 1) = 0, (3.4) 
u(0,x) — uq{x), (3-5) 

where we used the following notations 

ait) ^ 4-. m - ^^^^^'(^^;^r'^- (3-6) 
^' R^it)' ^' K2PGit,Rit))~K3 ^ ' 

Our purpose in this Section is to analyze the boundary value problem 
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We shall do so by setting this problem in a variational framework, using weighted Sobolev spaces. 
Let =]0, 1[ and define H as the Hilbert space given by 



H = {v.n 



measurable and such that / x'^v'^{x)dx < +00} 



Note that H is the closure of C^{n) w.r.t. the norm \\v\\h = 



x'^v'^{x)dx 



1/2 



We also introduce 



the real Hilbert space In the following, we shall often use the fact that 

„ „ /„ _ , \l/2 



V 



is the closure of C^{Vl) w.r.t. the norm \\v\\v = \ + . V is continuously embedded in 

H . Identifying H with his dual H\ one has V <Z H <Z y'with continuous injections. 
Note also that the norms and can be defined, respectively, from the inner products 



<u,v >— I x^u{x)v{x)dx and < u,v > + < u' ^v' >. 
We then have the following results, the proofs of which can be found in the papeA^ 
Lemma 3.1. For every v E C^{[0, 1]); e > and x E [0, 1] we have 



vHl)<e\\v'\\l + CM\l 



v{l) 



<2 Iklli, 



XV {x)\ < %/5|| 



where — S + t and 



Lemma 3.2. The embedding V d H is compact. 



Remark 3.1. From Lemma [3?T 



norms on 



it follows that ( ||w'||q + ti^(l) j and are two equivalent 



□ 



□ 



V since < ^^(l) + < 5||«||^, for aU veV. 



□ 



Remark 3.2. We have xv{x) S ^([0, 1]), for aU v eV. 

Indeed, on one hand, lim xv(x) — 0, ^v £ V (see the boolff, p. 128), and on the other hand 

e C°([e,l]), Ve, < e < 1, since we have H'^{e,l) C C°([e,l]) and e\\v\\Hi(,s) < Mi ^v e V 
Ve,'o<e<l. 

□ 

li X is any Banach space, we denote by ||.||x its norm, and by X' the dual space oi X. We denote by 
LP(0,T; X), I < p < 00, the standard Banach space of real functions u : (0,T) X, measurable, 
such that 
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\Lp((),T;X) 




1/p 



m'xdt 



< 



for 1 < p < oo 



and 



M\l--{o,t,x) =csssupo<t<Tll'«(t)|lx' for P = 00. 
Let u'(t), u^{t), u^^{t) denote u{t,x), ^{t,x), ^{t,^), respectively. 

We shall make the following set of assumptions: 

• (Hi) Mo e H; 

• (H2) a,ke W^i'°°(0, T), a{t) > oo > 0; 
. (Fi) /eC(M,R); 

• (F2) There exists positive constants Ci,C[, C2 and p, 1 < p < 3, such that 

(i) ufiu)>C,\u\P-C[, 

(ii) \fiu)\<C2il + \u\P-^). 

Let u e C2([0,T] X [0, 1]) be a solution of problem (1331) -(133]). 

Then, after multiplying equation (|3.3p by x^f , u G w.r.t. the scalar product of H , integrating by 
parts and taking into account boundary condition given by p.4p . we get 



— < u(t), t; > + a(t) / X UxVxdx + ait)kit)u{\)v{l) — / x Uxvdx+ < f(u), v >= 

at Jo R[t) Jo 

The weak formulation of the ibv problem (|3.3p - p.5p can then be given in the following way: Find 
u{t), defined on the open set (0,r), such that u{t) satisfies the following variational problem 

^ <u{t),v > + &{t;u{t),v)+ < f{u{t)),v >^0, yveV, (3.7) 
together with the initial condition 

u{0) = uo- (3.8) 
Above, we have used the following bilinear form 

i?'(t) 

a{t;u,v) = a{t) x^UxV^dx + a{t)k{t)u{l)v{l) / x^Uxvdx, u,v£V. (3.9) 

Jo Ja 

We first note the following lemma, the proof of which can be found in our previous papeJ^l 
Lemma 3.3. There exists constants Kt, otT and (3t depending on T , such that 

I a[t;u,v)\ < i^TllMllillwlli, for all u,veV, (3.10) 
d{t; u, u) > aT||u||i ~ I3t\\u\\o, u,v G V. (3.11) 
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□ 

We then have the following existence theorem 

Theorem 3.1. Let T > and assumptions (Hi),(H2),(Fi),(F2) hold true. Then, there exists a 
solution u of the variational problem {3. 7|), \3. such that 



u G 2.2(0, T; V) n L°°(0, T; H), x^l'^u € ^^(Qt), 

tu e L°°(0, T; V), tut e L2(0, T; H). 
Furthermore, if f satisfies the additional condition 

{f{u)^f{v)){u-v)>-S\u-v\', 
for all w, u G R, for some (5 G R, then the above solution u is unique. 



□ 



Proof of Theorem 13.11 We divide it in several steps. 

• Step 1, Galerkin method. 

Denote by {wj}, j = 1,2, an orthonormal basis of the separable Hilbert space V. We wish to 
find Um{t) of the form 

m 

Umit) ^^Cmj{t)Wj, (3.12) 

i=i 

where Cmj{t) satisfy the following system of nonlinear differential equations 

< u'„(<), Wj > + d{t; Umit), Wj)+ < f{u„i{t)), Wj >= 0, 1 <j < m, (3.13) 
together with the initial condition 

Um{0) = Uom, (3-14) 

and 

uom uo strongly in H. (3.15) 

Clearly, for each m, there exists an unique local solution Um{t) of the form p.l2p . which satisfies 
([Xrai and ([XTI)) almost everywhere on < i < T„, for some T^, < T™ < T. The following 
estimates allow us to take = T for all m. 

• Step 2, A priori estimates, 
(a) First estimate. 

Multiplying j*'' equation of system (j3.13p by c„ij{t) and summing up w.r.t. j, we have 

^^l|wmWllo+ ~a{t;u^{t),u^{t))+ < f{u^^{t)),Ura{t)>=Q. (3.16) 
Using assumption (H2), (F2,i), Lemma [3Tl and Remark l3.ll it follows from (|3.16p that 
d 2C 

-\\Ura{t)\\l + 2aT\\Um{t)\\l+2Ci j | U„ , x) l^dx < ^ + 2/3t || (i) || q • (3.17) 
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Integrating (|3.17p . using (|3.15p . it follows that 

2 
3 

where 



Sm{t) <Co + ItC[ + 2I3t I S^is)ds, (3.18) 







S^t) = \\um{t)\\l + 2aT f \\u^{s)\\lds + 2Ci f ds f x^\um{s,x)\Pdx, (3.19) 
Jo Jo Jo 

and Co is a constant depending only on uq with ||Mom|lo ^ Vm. 
Applying Gronwall's lemma, we obtain from (j3.18p 

S„^{t) < (^Co + ^TC;^exp(2/3Tt) <Mt, Vm, yt,0<t< T„, < T, (3.20) 
that is Tm = T. 

In the following, we denote by Mt any generic constant depending only on T. 

(b) Second estimate. 

Replacing wj by t'^Um in p.8p gives 



|, / ||2 l_rf 

' 2dt 



a{t)\\tu„^\\l + a{t)k{t)t^uUl) 



ld_ 
2dt 



• 1 

/ x'^f{u„i)dx 







= ll"m.|lo^ [tMt)] + \ul^{l)j^{t^a{t)k{t)) + 2tJ x'f(u„,)dx (3.21) 

-R(i) J ^rnxU^dx 

where 

f f{y)dy. (3.22) 
Jo 

Integrating (|3.2ip w.r.t. time variable from to t, we have, after some rearrangements 

2 / \\su'„^{s)\\lds + a{t)\\tu„,S)\\l + 1) 
Jo ^ 

= a{t){l~k{t))ful,{t,l)+ [ [s^a{s)]'\\u,nx{s)\\lds 

(3.23) 



+4/ sds / X f{um{s, x))dx ~ 2t / x f{um)dx. 
Jo Jo Jo 

By means of assumption (H2) and Remark l3.11 we get 

a{t)\\tu^{t)\\l + a{t)t^ul,{t,l) > ^aoWtu^ml Vt G [0,r], Vm. (3.24) 







+ [ [s^a{s)k{s)]'ul^{s,l)ds + 2 [ < xUr^^is) , u'^{s) > ds 



Jo 



We fix e > such that 



where ||.||oo = |M|loo(o,t)- 



klloolli-^lloo^<y. (3-25) 
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Using again Lemma 13.11 Remark 13.11 with e > as in l|3.25p and first estimate p.20p , the terms on 
the r.h.s. of (|3.23p can be estimated as foUows 



i(t)(l - k{t))fulit, 1) < ||a||^||l - k\\^[e\\tu,r.ml + CSumim 



< ^\\tUmit)\\l +Mt, 



a{s))'\\umxis)\\lds + / (s^a(s)fc(s))'u^(s, l)ds 



< [||(iV|L + ||(iW 



\\u,nx{s)\\Q+U^n{s,l) 
t 



ds 



^ 5[||(i'a)'|L + ll(*'«^)'IIJ / \Ms)\\lds < Mt, 

Jo 

2 „t 



m 



< XUmx{s),u'^{s) > ds 



< 



\su'^is)\\ods - 



ooJO 



From assumptions (Fi) and (F2), we note also that 



mo = 



\fiy)\dy < f{z) - f{y)dy < C2 {\z\ + ) ' ^ 



where zo = (C(/Ci)1/p. 

Using first estimate p.20p . (|3.29p and Lemma [3. 11 we obtain 

4/ sds I f {um{s , x))dx — 2t^ I x'^ f{um{t,x)dx 



< 4C2 / sds X \ \Urn,{s, x)\ -\ \Um(s,x)\ ]dx + 2t / XITlodx 

Jo^ Jo \ F 1 ^ ° 

< 4C2 / s\\um{s)\\r.ds + -C2t I ds x'^\umis,x)\''dx + -T^rho 

Jo P Jo Jo 3 

2C 2 
< 2C2TVMT + — -TMt + -T^mo < Mt- 
pCi 3 

Hence, we deduce from (IX^ . (021, and ((X^ that 



\\su'„,{s)\\lds + ^\\tuUt)\\l < Mt + 



2 



\sUm{s)\\ids. 



By Gronwall's lemma, we get 



\\su'„As)\\lds + y||tM„(i)||^ < Afrexp] 





R' 




R 



3T 
ao 



<Mt, yte[o,T]. 



(3.26) 



(3.27) 



lsu™(s)||ids. (3.28) 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



Finally, using (|3.20p and assumption (F2,ii) we have also 



ds [ 
Jo 



X^'P' f{um{s,x))^ dx<{2C2Y' ds x^\um{s,x)fdx<MT, 



with n' — — 

• Step 3, the limiting process. 



(3.33) 
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From (|3.20p . (|3.32p and (|3.33p . we deduce that there exists a subsequence of {um}, stiU denoted 
{um} such that 

u„,->u weakly * in L°°(0, T; H) , 
Um — > u weakly in L^{0, T; V) , 

^ weakly in Lf(QT) , „ 

tum tu weakly * in L°°(0, T; 1^) , ^ ■ 

(iwm)' -> (tu)'weakly in L2(0,T; ff). 

Using a standard compactness lemmsff^ {p-57) together with (|3.34p . we can extract from the sequence 
{um\, a subsequence still denoted by {m„i} such that 

tum -> tu strongly in L^(0, T; H). (3.35) 

Continuity of / also implies (up to a sub-sequence) 

f{um{t.x))^ f{u{t,x)) a.e. (t,a;) e Qt = (0,T) X (0,1). (3.36) 



Applying a standard weak convergence lemmsff^, we have also 

X^/P' f{um) X^/P',f{u) weakly in Lp'{Qt). 

Passing to the limit in ((XTO)) and ([XTI)) . it follows from (|XTS)) . (PIM)) and that function u{t) 

satisfies the i.b.v. problem p.7p . (|3.8p . 

• Step 4. Uniqueness of the solutions. 

First of alb we note the following slight extension of a lemma used in our previous papeJ^l (see also 
the boolJI^ 

Lemma 3.4. Let w be the weak solution of the following i.b.v. problem 
wt - a{t){wxx + x'^x) = f{t, x), <t <T, < X <l, 

lim xwx(t,x) < +00, Wx(t,l) + k{t)w{t,l) — 0, w{0,x)=0, 

x—*0+ 

w e L^{0,T-,V)n L^{0,T-H), x'^l^w e LP{Qt), 
tw e L°°iO, T; V), twt e L2(0, T; H). 

Then 



\\\w{t)\\l+ j a{s)\\\wx{s)\\l + k{s)w''{s,l)\ds~ [ < f{s),w{s) > ds = 0, 

a.e. t e (o,r) 

□ 

Uniqueness of solutions for our initial i.b.v problem will then be deduced as follows. Let u and v 
be two weak solutions of (|3.3p - (|3.5p . Then w — u — v is a, weak solution of problem mentioned in 

m-m ~ xR'(t) 

Lemma 13.41 with r.h.s. given by f{t, x) = p/)\ ~ /("") + /(w)- Therefore, Lemma 13.41 implies 

K[t) 
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^"w{t\\l+ [ h{s;w{s),w{s))ds + 2 I < f{u{s)) - f{v{s)),w{s) > ds ^ 0. 



2 

Using Lemma 13.31 and assumption (F3) we obtain 







\\w{tX + 2aT \\w{.s)\\i<2{d + l3T) \\w{s)\\ids. (3.37) 
Jo Jo 

If (5 + /3t > we have ||w(t)||o = by applying Gronwall's lemma. In the case where S + /3t < the 

result is clearly still true. 

This ends the proof of Theorem 13.11 

□ 

We now turn to the boundness of the above solutions. 

For this purpose, we shall make use of the following assumptions 

• (H'l) uo e i°°(0, 1), \uo{x)\ < M, a.e. x £ (0, 1) 

• (H'2) a,ke VKi'°°(0, T), a{t) > oo > 0, k{t) > fco > 

• (F'l) uf{u) > Vw e M such that |m| > ||moI|oo, for a.e. x S (0, 1). 
We then have the following result 

Theorem 3.2. Let (H'l), (H'2), (Fi)-(F^) and (F\) hold. Then the unique weak solution of the 
ibv problem ll3.7^ - [37^} . as given by theorem 1, belongs to L°°{Qt)- 

□ 

Proof of Theorem 13.21 Firstly, we note that Z ^ u — M satisfies the i.b.v. problem 

dtZ - ait) (az + -d^z] - x^^^d^Z + f(Z + M) = 0, < a; < 1, i G (0, T), (3.38) 
V ^ J 

\lim,^o+xZ,{t,x)\ Z,(t,l) + k{t)[Z{t,l) + M]^0, (3.39) 

Z{0,x) =uo{x) ~ M (3.40) 

Multiplying equation (j3.38p by x'^v, for v (E V, integrating by parts w.r.t. variable x and taking into 
account boundary condition (|3.39|) . one has 

x'^Ztvdx + a{t) [ x'^Z^v^dx + a{t)k{t)Z{t, l)w(l) - ^^7^ / x^Z.^vdx 



m 

+ I x^f{Z + M)vdx = -Ma{t)k{t)v{l), Vw £ V, 
Jo 

hence for v ^ Z+ = ^(^Z + since uq e L°°(0, 1). It follows that 
^-^J^ x^\Z+\^dx + a{t) x''\{Z+U^dx + a{t)k{t)\Z+{t, 

j ^^^xZ^dx + J x^f{Z+ + M)Z+dx = -Ma{t)k{t)Z+{t,l)<0, 

since 

f x'ZtZ+dx^ C x\Z+)tZ+dx^~ f x^\Z+\^dx. 
Jo Jo.z>o 2 dt Jq 



(3.41) 



15 



On the other hand, by assumption (H'2) and Remark 1 3. 1[ one has 



(t) x^\Z+\^dx + a{t)k{t)\Z+{t,l)\^ > Co\\Z+{t)\\l, 



(3.42) 



where Co — min {1, fep}. 

Using the monotonicity of f{u) + Su and (F'l), we have 



x^f{Z+ + M)Z+dx = / x^[f{Z+ + M)-f{M)]Z+dx+ / f{M)x^Z+dx 
^ Jo ^ Jo 

>-s[ x^\Z+\^dx+ [ f{M)x^Z+dx > -S\\Z+\\l. 
Jo Jo 



(3.43) 



p.4H) - (|3.43p together with Cauchy's inequahty apphed to the term 



m Jo 



-\\z^{t)\\;+Co\\z^t)\\i < [- 



Co 



+ 2\S\]\\Z+{t)\ 



x^Z+Z+dx yields 



(3.44) 



Integrating (|3.44p . we get 

||Z+(t)||^<||Z+(0)||J + 



Co 



+ 2\5\\ \\Z+{s)\\lds. 



(3.45) 



A/ ) ^ = 



(^0(2;) - ^ = 0, Gronwall's lemma yields ||^+(<)||o = 0. Thus 



Since Z+(0) = {u{Q,x) 
u{t,x) < M a.e. {t,x) £ Qt- 
The case uo{x) > —M is similar, by considering Z = u + M and Z~ = ^(\Z\ — Z^ . Thus we get 
Z~ = and hence u{t, x) > —M a.e. (t, x) G Qt- 

All in all, one obtains \u{t, x)\ < M a.e. {t, x) € Qt and this ends the proof of Theorem 13.21 

□ 



4. Numerical applications 

For the numerical applications, we have taken in (|2.3p 7 = 3, so that equation (12. 9p reduces to 
Burger's equation 

(Wt + WWr^O, W{0,r)^Wo{r) 

\ Zt + ZZr^O, Z(0,r) = Zo(r). ^ ' ' 

It is well known that classical Burger's equation 

Ut + uUr — 0, u(0, r) = UQ{r) 

admits the solution u{t, r) — uo{S,(t, r)), ^(i, r) being defined by the parametrization r = uo{^)t + 
Having in mind (|4.ip . we have considered two examples. 
• First example. 

For the first example, we have chosen the initial conditions Wo{r) = 1, r > 0; Wo(0) = and 
Zo{r) = 2, r > 0; Zo{0) = 0. The continuous solutions of (|4.ip are then given by 
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W{t,r) 



1 if < t < r, 
f if < r < 



2 if < 2t < r, 
? if < r < 2t. 



According to Section 2, the droplet radius is given by formula (|2.19p which we consider here with 
an initial condition taken equal to be 1 



dR{t) VG{t,R{t))pG{t,R{t)) 



dt 



and where 



PG{t,R{t))- PL 



1 



i?(0) = 1, 



VG{t,r) = -{W{t,T) + Z{t,r)) 
PG{t,r) ^ ^^^{Z{t,r) -W{t,r)). 



(4.2) 



(4.3) 



In figure 1 below, we have drawn the curve t — > R{t) on the time interval [0, 1] with a step h ~ 0.05 
and PL = 0.9. 




Fig. 1 

• Second example. 

For the second example, we have chosen the truly experimental conditions made by the LCSR in 
the study of single drop evaporation, the drop being suspended from a silicate tube. Drops are 
made up of n-heptane fuel {pL — 683 kg/mm?)) in air at normalized atmospheric pressure and 
with an initial speed wg(0, r) = Ci — 35mm/s. The initial density pG(0,r) of the gas is taken as 
C2 = To = 373K. 

In this case, the solution of (|4.ip are given by 



1 + J 1 - A{r - Cit)V3C2t 



2%/3C2t 
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1 + Vl + 4(r-Cii)\/3C2i 



We then compute vg and pc by formula (j4.3p . and then solve the ode for R given by 



We note that R'{0) = < where C2 ^ ^ < K Pl = 683. 

On the other hand, one has 



voit, r) = ^[C, + + Vl—a - ./TTT^)], 



where 



a = a(t, r) = 4(r - Cit) V3C2i, |a| < 1. 



One has sgnwcl^,^) — sgn[C\12t^C2 + V3(a + \/l — a — \/l + a)]. Since |q;| < 1, it follows that 
V3(a + Vl - a - %/l + «) > -Ve. 

Thus WG(i, r) > if Cil2i2C2 > V6 i.e. t > t2 = ^-r =i ^ and vcit, r) < on the interval [ti, • 

Similarly < pait^r) < dA ^.r^ [2 + VI - a - VI + «] < 34^- 

Thus, if t > t„j, r > r„i, pait^r) < 3^/^^ < Pl ^ 683, that is for i„r„ > ~ ^. 

Since i?'(0) < it follows that R{t) is decreasing on [0, increasing on [ti,t2[ and then from the 

starting point ^2 always non increasing. 

Figures 2 and 3 represent resp. the velocity vcit^r) and the pressure pG{t,r) given by (4.3) for 
(i,r)g (0,l)x (0,1). 




Fig. 2 
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Fig. 5 

Since Wq (^) < the maximal existence interval is finite (Prop ositi on 12. ip as can be seen in our 
graphic. Let us remark that looking on the experimental curves made by the LCSR (figure 5) 
at the beginning, the function t — > R{t) is increasing around the vicinity of t = 0. This fact is 
confirmed by our model which represents a good improvement of our previous model^ in which the 
velocity VQ{t) was a given function of t. 
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